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Abstract 

The satellite SEASAT-A will carry a radar scatterometer 
in order to measure microwave backscatter from the sea 
surface. From pairs of radar measurements at angles 
separated by 90° in azimuth the surface wind speed and 
direction may be inferred, though not uniquely. In this 
paper the character of the solutions for wind speed and 
direction is displayed, as well as the nature of the ambi- 
guities of these solutions. An economical procedure for 
handling such data is described, plus a criterion for the 
need for conventional (surface) data in order to resolve 
the ambiguities of solutions. 
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THE NATURE OF MULTIPLE SOLUTIONS FOR 
SURFACE WIND SPEED OVER THE OCEANS 
FROM SCATTEROMETER MEASUREMENTS 

I. INTRODUCTION 

A knowledge of surface wind speed and direction over the oceans is desireable 
on a number of counts. The value to shipping and fishing industries and to the 
military is quite apparent. Predictions of winds and waves are likewise im- 
portant for many eoastal applications, with the case of a hurricane making 
landfall an obvious example. 

In a broader sense surface wind is an ingredient in the understanding of the 
global circulation. Through surface friction the wind provides the driving 
force for small and large scale ocean currents: it is thus an essential input 
to ocean dynamics models. For meteorological applications the knowledge of 
the spatial distribution of surface winds provides information on the surface 
pressure distribution, which in turn is a key parameter ‘V « global weather 
forecasting. In all these cases the need beyond current . <, « ailable data is for 
global measurements at frequent time intervals. 

The radar scatterometer to be flown on SEASAT-A has generated considerable 
interest because it provides data which may be analyzed to yield the desired 
winds. Individual measurements, having a footprint of some hundreds of square 
kilometers, yield estimates of spatially averaged wind values, with the local 


(point) variability due to guatineaa averaged out. These values are the required 
input for oceanic and atmospheric applications studies. 

The purpose of this paper is to describe the interpretation of radar data, and 
to develop an efficient procedure £« v obtaining values for surface wind speed 
and direction from the measurements. 

The general features of the measurement system and of the nature of the radar 
return are described in section II, which draws heavily on work by Pierson, 
Cardona and Greenwood (1974). Section Hi describes a lookup table procedure 
which may be used in "inverting'* the model equer -ons for radar backseatter 

to obtain wind speed and direction, as well as the nature of the solutions. In 
most cases there are four possible values for the wind angle, which may be 
expressed as either a, 3(j<j u -m and 180° ±0 , or a, 27(T-a , and SU' ±0, with a, 
0< 90*. The necessity for conventional data to resolve this ambiguity is 
discussed in section IV. 

H. S C ATTGROM E TER MEASUREMENTS 

Although the SEASAT is not completely defined at this time, a general descrip- 
tion of the seatterometer is possible. From a height of approximately 800 
kilometers the instrument will measure radar backseatter at 13.9 GHfc ( X 35 2.2 
cm). Measurements will be taken in a crosstrack pattern, with radar foot- 
prints at singles (in azimuth) 45° to right and left, fore and aft of the subsatel- 
lite track. Due to satellite motion the backseatter from a spot on the ocean 
surface will be measured twice, at nearly coincident times, but with a 



separation of 90° In the direction of measurement (see Figure 1). These two 
measurements of radar cross section o , permit estimation of two quantities, 
surface wind speed V, and direction x • This geometry appears to be optimum 
for observation from a satellite, given the level of instrumentation which is 
currently feasible. Observations at a number of azimuth angles would be 
deaireable, but this would cause greater complexity in the radar scanning 
mechanism and electronics. The dependence of radar cross section on surface 
roughness lias been estimated theoretically (Beckman mid Spizaaehino (1903), 
Stogryn (1967), Wu mid Fung (1972 i) but experimental results are not completely 
understood. For present purposes the dependence of cross section on surface 
wind is assumed, and the problem is that of inverting c^, (x . V), o 2 ( % + 
9U°, V) to infer x and V. 

By scanning in nudir angle the seatterometer will provide wind field estimates 
in swaths on both sides of the subsatellite track. (At near nadir viewing the 
pairs of measurements are not independent and no information is available on 
wind direction). Questions of instrument calibration, noise, etc. are not 
addressed here as these have no direct bearing on the development of mi ef- 
ficient procedure for converting measured radar cross sections into surface 
winds. 


A starting point is the equation for radar backscatter derived by Pierson et al 
from experimental data from the Langley A A FE (Advanced Applications Flight 
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Experiment) Program. For vertically polarized transmitted and received 


radiation this equation; is 

<T VV ( x , V)«* * (EV n - K 2 V^) COS2X + 1/4 (KgV ,U + K,V q - 2EV* 1 ) x 

(3 cos \ + cos 3 X i (1) 

where (E, K,, Kj, N, M, Q) are constants depending only on the nadir angle, 

V is the wind speed, and x is the angle measured clockwise from the meteoro- 


logical wind to the measurement angle. (The meteorological wind direction 
is that from which the wind blows). Pierson et al give values for E “-Q based 
on AAFE Kadseat data and Skylab data, both at 13.9 GHz. Preliminary values 


for a nadir angle of 30° are: 


E - 2, 24 x liT 8 

M = 1. 83 

K, - 2.27 ?: 10~ 3 

N- 1.69 

K 2 - 2. 21 x 10 -3 

Q - 1.55 


These values were used in this study in order to model the radar backsoatter. 
Revised values, or even a different functional form (e. g. , including terms in 
cos4 X ) may be necessary as more and better data are analyzed. Such changes 
will not affect conclusions drawn here ay long as the general behavior of the 
radar cross section is not changed. Figure 2 illustrates the dependence of 
a on x and V, with values given in decibels due to the great dynamic range 
involved. Henceforth a will be referred to in db. 
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irt. SOLUTIONS IT OR WIND SPEED AND ANGLE 


The proposed sut6''v,e instrument will obtain nearly simultaneous measure- 
ments of cross section u at angles separated in azimuth by 90°, i. e. , % 

(X , V), (\ + 90°, V). Although the solution for x , V of the equations 


*i C -V , V) * 


, measured 
a i 


a 2 ( X + 90°, V) « 


measured 

a 2 


is reasonably straightforward, fox* large scale data reduction or for operational 
purposes a lookup table is preferable, From tills table questions of the am- 
biguities of solutions may be worked out systematically, and computer savings 
are possible once the lookup table is generated. 


In this study it was found that the quantities ° l and - °i may be used to label 
a storage array with a minimal number of entries which includes the full range 
of values of c 2 to be expected. The first of these, Oj, provides a reason- 
able first estimate of the wind speed , while the difference, u 2 ~ °i» provides 
information principally about the wind angle . Because the radar backseatter 
tends to increase faster with windspeed at X - U°, X - 180° than at X - 9(J U , 

X - 270° (see Figure 2) the magnitude of a 2 ~ ^ increases with surface wind- 
speed, and thus with ^ For this reason it is best to scale » 2 ~ °i by dividing 
by 

D= 4. 50 + 0.084 ( °i + ff 2 ), 
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which is to good approximation the minimum value of Oj. Uj, for a given value 
of 0 |, It follows that the ratio 5 * (u 8 - 'ql/D varies between ~ 1 and 1, and 
tables of V (o t , 6 >, X ( </ it 6 ) represent all physically possible combinations 
of ffg in a compact fashion, granted t e hypothesis that equation 3 adequately 
describes the measurements. 

The non uniqueness of X and V is Illustrated in Figures 3 through G. Values of 
have boon chosen in order to represent regimes of high wind speed ( ‘J t «* 
-3.5 db) find low wind speed («q « - 10. U db). The degree of ambiguity of a 
solution is giv&n by the number of times a vertical line ( fi - constant) inter- 
sects the solution curves. The solution regimes are most easily described in 
terms of the angle x between the first measured value aud the meteorological 
wind. 

1. Values of S from - 1. 0 to approximately - 0. 7 represent a "first 
look upwind" solution in which the first measurement is taken 
nearly upwind, the second nearly cross wind. The results illus- 
trate the inability to distinguish left-right values of angle, i. e. , 
u , 3G0°- w, with 0<a<45°. 

2. Values of 8 from i. 0 to approximately 0. 7 represent a "second 
look upwind" solution in which the first measurement is approxi- 
mately cross wind, the second is nearly upwind. The resultant 
ambiguity with respect to left -right of upwind corresponds to 270° 
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± u us defined by the angle of the first radar measurement with 
respeet to the wind. 

3. For all intermediate values of 6 7 to+ , V) there exist, in addi- 

tion to "upwind" solution 8 ut larger angles relative to the wind 
direction, the "downwind" solutions which begin at approximately 
7 - "first look downwind, " and + 0. 7 - "second look downwind. " 
These values of 3 correspond to angles from 45° to 2US U , 

The results of ambiguity in wind speed determination are easier to describe. 
Basically the solutions fall into pairs representing upwind and downwind. As 
illustrated in Figures 5 and V the upwind solutions fall on the lower, larger 
wing of the butterfly shape, and the upwind solutions fall on the upper, smaller 
wing. 

It should be emphasized at this point that the parameter 8 involves the dif- 
ference in the logarithms of the measured radar returns, and is thus sensitive 
to small values due to the derivative d (log x) = dx/x. Small values occur 
when either measurement is crosswind. Such eases arc particularly bad, 
because such values also represent turning points of the solution for wind angle 

oo) j n 3 and 4, in such a ease some form of data averaging 

d b 

or smoothing is absolutely required. The next section recommends such a 
procedure. 
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IV. RECOMMENDATIONS FOR DATA ANALYSIS 


The previous section illustrates the fact that seatterometer data, by ite?,!*, 
does not yield satisfactory (unique) results for surface winds over the oceans. 
Additional input data will be required* such as surface wind observations, or 
results from a global scale meteorological analysis, which relies on observa- 
tions of surface wind and barometric pressure, in either case the spatial 
density of actual surface measurements will be much lower than the density 
of seatterometer measurements. For physically reasonable results it will be 
necessary to assume spatial continuity of wind fields on tile settle of 2U-5U 
V=' - neters, and to average or smooth the scaiterometer results. The assump- 
tion of spatial smoothness is not a strong one, as synoptic scale wind variations 
are generally significant only over distances of many hundreds of kilometers. 

In principle one may smooth either the measurements Uj, u a of the seatter- 
ometer, or the resultant values of derived wind. However the sensitivity of 
results for wind angle at o°, 90°, 180°, 270° to variations in measured cross 
section argues strongly in favor of smoothing the seatterometer data itself. 

The degree of such smoothing will represent a trade off between effects, in- 
strument noise errors and the amount of spatial smoothing which is acceptable. 
Undoubtedly experimentation on the amount of smoothing or averaging will 
have to be carried out after the launch of the satellite. It is possible that 
different averaging routines will be needed for different purposes. 
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From the analysis it is possible to specify the volume of ground truth required 
in order to specify unique wind fields from the scatterometer results. In an 
approximate sense the solution for wind angle is given by \ ~ -cos -1 (2 * ). 
From the spatial continuity of surface winds it follows that an ambiguity must 
be resolved only in the vicinity of possible movement from one branch of the 
arc cos function to another, i. o. , at u°, 90°, 1SU° and 270°. In terms of an 
automated analysis of scatterometer data, only these branch points must be 
flagged. After such identification ail intermediate data may be specified by 
comparison with surface truth data in the same region. As little as one surface 
measurement may be used to fill in uniquely the wind values from one branch 
point to the next. 


Such a procedure is considerably simpler than that advocated by Pierson, 

Cardona and Greenwood, which involves comparison with all available high 
quality ship reports, with appropriate interpolation to points midway between 
ship reports. Such a procedure is excessively conservative, provided dy 
that the scatterometer data is of adequate quality to permit reasonable smooth- 
ing. Of course the analysis described here of comparison with ship reports 
once per solution branch does not preclude further Comparison with all avail- 
able surface observations. Hopefully such comparisons will soon establish 
the validity of the satellite derived values, and the resultant global wind field 
over the oceans. At this point a compact lookup table plus occasional reference 
to surface truth will permit rapid and efficient application of the SEASAT results. 
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Figure 1. View From Above of Scatterometer Measurement Geometry. The In- 
strument is Planned to View Left as Well as Right (Shown), and to Scan in Nadir 
Angle from 25 to 45 , Providing Parallel Swaths Along the Suhsatellite Track 



Figure 2. Variation of Radar Crons Section Versus Look Angle for Represent- 
ative Values of Wind Speed. The Curves are Symmetric About 0 
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Figure 3. Solution Curve for Wind Angle Moderately High Winds ( = -2. 5) 

as a Function of the Difference in Radar Cross Sections 



Figure 4. Solution Curve for Wind Angle for Low Wind Speed (o, = -10. 0) as a 
Function of the Difference in Radar Cross Sections 
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(«2 * ®i)/D 

Figure 5. Solution Curve for Wind Lpeed for High Winds as a Function of the 

Difference in Radar Cross Sections 
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Figure 6. Solution Curve for Wind Speed at Low Winds (<*, - -10. 0 db) 
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